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Abstract 

As a main result of the paper, we construct and justify an asymptotic approximation of Green's 
function in a domain with many small inclusions. Periodicity of the array of inclusions is not required. We 
start with an analysis of the Dirichlet problem for the Laplacian in such a domain to illustrate a method 
of meso scale asymptotic approximations for solutions of boundary value problems in multiply perforated 
domains. The asymptotic formula obtained involves a linear combination of solutions to certain model 
problems whose coefficients satisfy a linear algebraic system. The solvability of this system is proved 
under weak geometrical assumptions, and both uniform and energy estimates for the remainder term are 
derived. 

In the second part of the paper, the method is applied to derive an asymptotic representation of 
the Green's function in the same perforated domain. The important feature is the uniformity of the 
remainder estimate with respect to the independent variables. 

Keywords: Singular perturbations, meso scale approximations, multiply perforated domains. Green's function 

1 Introduction 

Uniform asymptotic approximations of Green's kernels for various singularly and regularly perturbed domains 
were constructed in [l]-[4]. In particular, the papers [1], [2] address the case of domains containing several 
small inclusions with different types of boundary conditions. In the present paper, a similar geometrical 
configuration is considered, but the number of inclusions becomes a large parameter, which makes asymptotic 
formulae in [1], [2] inapplicable. 

In Sections 4-7, we address the Dirichlet problem for the Poisson equation — Au = / in a multiply 
perforated domain with zero Dirichlet data on the boundary. Section 4 contains the formal asymptotic 
representation 

N 

u(x) ~ i;/(x) + ^C'j(p(^'Hx) -47r cap(F(^)) i/(x, O^^'^)) , (1) 
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where 



• Vf is the solution of the same equation in a domain O without inclusions, 

• P^^^ is the harmonic capacitary potential of the inclusion F^^^ , 

• cap(-F(''^) is the harmonic capacity of F^^\ 

• H is the regular part of Green's function G of CI. 

The coefficients Cj should be found from the algebraic system 

(I + SD)C + V/ = 0, (2) 
where I is the identity matrix, and the matrices S,D and the vectors C, V/ are defined by 

S = {(l-(5,fe)G(0W,0«)} , D = 47rdiag{cap(F(i)),...,cap(FW)}, (3) 

and 

C = (Ci,...,C^f, yf = {vf{0^'^),...,Vf{0^^y)f, (4) 

with O'-'^ being interior points in F^^h 

The unique solvability of the system (2) is not obvious, and it is established in Section 5 under the 

natural assumption e < c d, where c is a sufficiently small absolute constant, e is the maximiim of diameters 
of F'^^\j = 1,. . . ,N, and d characterizes the distance between inclusions. In the same section, we obtain 
auxiliary estimates for the vector C under two diff'erent constraints on e and d. 

Justification of the formal asymptotic approximation (1) is given in Sections 6 and 7. where we show that 
the remainder admits the uniform estimate 0{s + £^(i~^/^) and the energy estimate 0{s^d^'^). 

Although we see that the asymptotic method described above may be applied formally under a very mild 
geometrical constraint e < c d, the convergence of the approximation in the space Lqc (CI) and in the Sobolev 
space H^{fl) has been proved when e <€. d"^^^ and £ <C d^, respectively. Hence, the asymptotic approximation 
(1) is efficient for certain meso scale geometries, intermediate between a collection of inclusions whose size e 
is comparable with d and the classical situation with e ~ const d^ appearing in the homogenization theory 
(see [5], [6] et al.). As is well known, in the latter case u is approximated by a solution of the equation with 
a "strange term" —Am + fj.u = f, jU > 0. 

In the concluding Section 8, we derive the above mentioned asymptotic formula for Green's function 
GAr(x, y), uniform with respect to x and y. The following is a specially simple form in the case of f2 = ffi.^ : 

^ - N ^ 

^^^'''^^ = 4dx-vl + Ef^'H'^'y) + E C,,J'W(x)P«(y) + 0{ed-'), 

where g^^^ are Green's functions in M.^\F^^\ and the matrix C = (Cij)^-^i is defined by 
C = (I + SD)-iS. 
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2 Main notations and formulation of the problem in the perfo- 
rated region 

Let fl be an arbitrary domain in R^, and let {O^^^jLi and {F^^^}jL-^^ be collections of points and disjoint 

compact subsets of fl such that O'^-') G and _F(^^ have positive harmonic capacity. Assume that the 

diameter ej of F^^^ is small compared to the diameter of fi. We shall also use the notations 



d = 2 ^ min 

ijtj,i<i,j<N ' 



max £ 

l<j<JV 
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It is assumed that e < c d, with c being a sufficiently small constant. 
We require that there exists an open set w such that 



N 

y F^^^ c oj, diam(a;) = 1, dist {duj,dn) > 2d, 
Let us introduce the complimentary domain 

as shown in Fig. 1. 



JV 

and dist{ |J F^^\doj^ > 2d. 




(5) 



(6) 



(7) 



Figure 1; Perforated domain containing many holes. 

Let u denote the variational solution of the Dirichlet problem 

-Aw(x) =/(x), xGfijv, (8) 
u(x) =0, X e dClN, (9) 



3 



where / is assumed to be a smooth function with a compact support in fl, such that diam(supp f)<C with 
C being an absolute constant. 

We seek an asymptotic approximation of u as A'' — > oo. 

3 Auxiliary problems 

We collect here solutions of some boundary value problems to be used in the asymptotic approximation of 
u. 

3.1 Solution of the unperturbed problem 

By vf we mean the variational solution of the Dirichlct problem 

-Avf{x) =/(x), xefl, (10) 
Vf(x) =0, xedn, (11) 

where / is the same smooth function as in (8). 

3.2 Capacitary potentials of F^^^ 

The harmonic capacitary potential of F^^^ will be denoted by P^^^ , and it is defined as a unique variational 
solution of the Dirichlet problem 

ApO')(x) =0 on R^\F^'\ (12) 
P(^)(x) =1 for X e a(R^ \P(^'), (13) 
P(^')(x) =0(£|x-0(j')|-i) as e-i|x-0(^')| -^oo. (14) 

It is well known (see, for example, [7]), that these functions have the following asymptotic representations: 
P(^')(x) = I'^^^Qtj)^! + 0{s cap(F(^'))|x - O'^'^r^) for |x - O'^'^l > 2e. (15) 
The harmonic capacity of the set F^^^ can be found by 

cap(F«) = -l/ |VP«(Opde (16) 

471" Jr3\fU) 

3.3 Green's function for the unperturbed domain 

Green's function for the unperturbed domain is denoted by G(x, y), and it satisfies the boundary value 
problem 

A,G(x,y)+<5(x-y)=0, x,yer!, (17) 
G(x,y)=0 as x e 917 and yen. (18) 

The regular part of Green's function is defined by 

iJ(x,y) = (47r|x-y|)-i-G(x,y). (19) 
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4 Formal asymptotic algorithm 

Let the solution u of (8), (9) be written as 

m(x) =t;/(x) + i?W(x), (20) 

where Vf solves the auxiliary Dirichlet problem (10), (11) in the unperturbed domain, whereas the function 
R^^^ is harmonic in Ojv and satisfies the boundary conditions 

i?W(x) = when X G 00, (21) 

and 

i?(i)(x) = -w/(x) = -t;/(0(*^)) + 0(£) when x e ^(IR^ \ F^'^)). (22) 
Let us approximate the function R^^^ in the form 

AT 

i?(i)(x) ~ ^C, (p(^)(x) -477 cap(F«) //(x,0(^))), (23) 
i=i 

where Cj are unknown constant coefficients, and P^-'^ and H are the same as in (12)-(15) and (19), respec- 
tively. 

By (15), (19) and (18), we deduce 

P(j^(x) - 47r cap(p(-''') i?(x, O^^^) = 0{e cap(p(-''')|x - O^^^T^), (24) 

for allx G dfl, j = 1,...,N. 

On the boundary of a small inclusion F^''^ [k = 1,. . . ,N) we have 

^;/(OW) + 0(£) + Cfc(l + 0(£)) (25) 

+ ^^"(^Tr cap(p(J')) G(0('=\0(J'^) + 0(£ cap(p(j'))ix - O^^'^i-^)^ = o, 

for all X G 9(R3 

Equation (25) suggests that the constant coefficients Cj, j = 1, . . . ,N, should be chosen to satisfy the 
system of linear algebraic equations 

^^(OW) + Cfc + 47r cap(P^^'^) G{0^^\0^''^) = 0, (26) 

l<j<JV, ji^k 

where k = 1 , . . . , TV. 

Then within certain constraints on the small parameters e and d (see (5)), it will be shown in the sequel 
that the above system of algebraic equations is solvable and that the harmonic function 

N 

P(2)(x) = p(i)(x) - ^C,(p(j')(x) - 47r cap(p(j')) i7(x, O^^'))) 

is small on 90 at. Further application of the maximum principle for harmonic functions leads to an estimate 
of the remainder P^^^ in Ojv- 
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Hence, the solution (20) takes the form 

N 

u(x) = i./(x) + ^q(p(J''(x) -47r cap(p(^)) H{x,0^^'>}'^ +i?(2)(x), (27) 
where Cj are obtained from the algebraic system (26). 

5 Algebraic system 

In this section wo analyse the solvability of the system (26), and subject to certain constraints on e and d, 
derive estimates for the coefficients Cj, j = 1, . . . ,N. 
The following matrices S and D will be used here: 

S = {(l-(5ife)G(OW,OW)} , (28) 

and 

D = 47r diag {cap(F(i)), . . . , cap(p(^^)}. (29) 

If the matrix I + SD is non-degenerate, then the components of the column vector C = (Ci, . . . , Cjv)^ are 
defined by 

C = -(I + SD)-iV/, (30) 

where 

V/ = (t;/(0«),...,t;/(OW))^. (31) 

Prior to the formulation of the result on the uniform asymptotic approximation of the solution to problem 
(8)- (9), we formulate and prove auxiliary statements incorporating the invertibility of the matrix I + SD 
and estimates for components of the vector (30). 

Lemma 1. // maxi<j<7v cap(p(-'^) < 5rf/(247r), then the matrix 1 + SD is invertible and the column 
vector C in (30) satisfies the estimate 

^cap(F«) < (1 - — ^max^cap(FW))-2^cap(FW) {vf{0^^^)f. (32) 

Proof: 

According to (30), we have (I + SD)C — — Vj. Hence 

(C, DC) + (SDC, DC) = -(V/, DC). (33) 
Obviously, the right-hand side in (33) does not exceed 

(C,DC)i/2(V/^DV/)V2. (34) 

Consider the second term in the left-hand side of (33). Using the mean value theorem for harmonic functions 
we deduce 

(SDC,DC) = (47r)2 ^ G(0(*\ 0(^')cap(p(*))cap(p(^') QC, 

i'?^j,l<i,j<N 

i^j,ikj<N mm Jb^Jbu^ 
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where B*^^^ = {x : |x — 0'^-')| < d}, j = \, . . . ,N, are non-overlapping balls of radius d with the centers at 
O^-'', and \B'^^^\ = 4nd'^/3 are the volumes of the balls. Also, the notation Bd is used here for the ball of 
radius d with the center at the origin. 

Let S(x) be a piecewise function defined on fl as 



S(x) 

Then 



CjCap(FW) in B'^i\j = 1,...,N, 
otherwise. 



(SDCDC 



V(cap(F«))2c2 / / G(X,Y)dXc^Yy (35) 

The first term in the right-hand side of (35) is non-negative, which follows from the relation 

G{X,Y)E{X)E{Y)dXdY = [ Vx / G(X, Y)S(Y)dY ^dX > 0. (36) 

Jn Jn 



The integral 

/ / G(X, Y)dXdY 
in the right-hand side of (35) allows for the estimate 

rfXfiY 



/ / G(X,Y)rfXdY< / / 
^^Ib, '^^^ /ykixi ^ i>|Y|>|x[ P^^} 

471" l-^O y{Y:|Y|=p} |X - Y| 7|x| y{Y:|Y|=p} |X - Y | J 

Using the mean value theorem for harmonic functions we deduce 

/ ^_=4V|X|-i when|X|>p. 



/{Y:|Y|=p} 

On the other hand. 



'{Y:|Y|=p} |X- Y| 

which follows from the relation 

dSy r d 



r dSy f 

Ay:IYI=„>P|X-Y| Jr 



■dSy 



{Y:|Y|=p} P|X - Y| y{Y:|Y|=p} d\Y\ |X - Y 

■/ Ay ^ dY = An when |X| < p. 

i{Y:|Y|<p} |X-Y| 



(37) 



(38) 



/ ,J\r> = 47rp when |X| < p, (39) 
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It follows from (37), (38) and (39) that 



/ / G(X,Y)dXdY< 1 / fd^-^^)dX 
Next, (33), (34), (35) and (40) lead to 
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(40) 



N 

(SDC,DC) >^jj G(X,Y)S(X)S(Y)dXdY - ^ ^^(cap(F(^")))2c|, (41) 



where a = ff • Then (33) and (41) imply 

9a 



15 • 

TV 



(l - ^ ^max^cap(J^(^))) ^C|cap(F(^)) 

(N \ 1/2 / AT \ 1/2 

^ C|cap(F(^)) j [^^(^KO(^)))^cap(F(^)) j , 



which yields 

/ N \ 1/2 / AT \ 1/2 

(l-^^max^cap(FO)))(^^C|cap(F«)J < (^5^(z;/(O«))2cap(F0)) j . (42) 

Thus, if maxi<j<Ar cap(F(-'^) < -^d, then the matrix 1 + SD is invertible and the estimate (32) holds. The 
proof is complete. □ 

Replacement of the inequality e < cd by the stronger constraint e < ccP leads to the statement 
Lemma 2. Let the small parameters e and d, defined in (5), satisfy 

e<cd^, (43) 

where c is a sufficiently small absolute constant. Then the components Cj of vector C in (30) allow for the 
estimate 

\Ck\ < c max \vf(0^^h\. (44) 



Proof: Let us write the system (26) as 

cap(f (j)) 

I rO) I 



Cfe+47r ^^^-^fSyr /o) G(OW,y)dy = -^;KOW), fc = l,...,7V, (45) 



where -Bj^^ is the ball of radius d/A with the centre at O^-'^ Also let cr be a piece- wise constant function 
such that 
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Multiplying (45) by cap(F('^' ) and writing the equations obtained in terms of a we get 

C7(0W) + — cap(i^W) / a(y)G(OW,y)dy = -Vf{0^''^) cap(F('=)), 



which is equivalent to 



a(OW) + i^ cap(i^W) / G(y,z)a(y)dy = cap(FW) $W(z), fc = l,...,iV, (47) 



where 



1 no /■ 

$('=)(z) = -VfiO^'^) + ^ y G(y,z)a(y)dy (48) 
/ a(y)((i/(0(*=),y)-i/(z,y))dy 



192 

48 



Next, we multiply (47) by 

/ N 2M-1 



G(y,z)a(y)dy 



<i/4 



where M is a positive integer number. Also, taking into account that a{0^''^) = a{z) for all z G we 
write 

2M-1 / X 2M 



/ G(y,z)a(y)rfy <7(z) + ^ cap(FW) / G(y,z)a(y)rfy 

(N 2M-1 
/ G(y,z)a(y)dy) , ^ & b'^^/,. 



Since a = outside the balls it follows that the integration of the above equation over B^^J^ and 

summation with respect to k = 1,. . . ,N lead to 

(X 2M-1 
J^G{y,z)a{y)dy\ a{z)dz 

N / \ 

= ^cap(i^«) /^^^$W(z) /G(y,z)a(y)dy dz. 



(49) 
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The identity 



2M-1 



^ I ^ G(y, z)a(y)dy j a{z)dz = 



M 



shows that the first term in the left-hand side of (49) is non-negative. By Holder's inequality, the right-hand 
side of (49) does not exceed 

/ N \ 1/(2M) 

5:cap(FW) /^^^($W(z))2MdJ 
\ fe=i ^^d/i J 

/ M ( N 2M X (2M-1)/(2M) 

x|^|^cap(i^W) 1^^^^ |^|^G(y,z)a(y)dyj rfzj 



and hence (49) yields 

2M N 1/(2M) 



192 



After the limit passage as M ^ oo we arrive at 



d ■ sup 
and by (47) we deduce 



G(y, z)a(y)rfy 



< c max sup |$W(z)|, 



(50) 



AT / \ \ . 

^ cap(FW) / / G(y,z)a(y)dy rfz 

.V \ 1/(2M) 

^cap(FW) / ($W(z))2^dz . (51) 



< c cap (i^^*^)) max sup |$(^)(z)|. (52) 

In turn, it follows from the definition (48) of the functions ^'^^^ that 

sup |$W(z)|<|«/(0W)| + i^ max KOW)| sup [ G{y,z)dy 

+ max |a(OW)| sup ^ [ \H{0^'^\y) - H{z,y)\dy 

+ ii max |a(0(9))| sup V / ^ ^^^^tSt^ jdy, 

.rf3 i<,<^l \e4,<,4:,^,iB<^, |y-z||OW -y| 

which, together with (52), yields 

|a(OW)| <ccap(FW) { m^ |^;/(0«)|+d-2 max^l^^^^^^ 

L 1<J<N 1<J<N ) 

If maxi<fe<;v cap(F('^)) < cd^, with c being a sufficiently small constant, then referring to the definition (46) 
of the function a we deduce (44), which completes the proof. □ 
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6 Meso scale uniform approximation of u 

We obtain the next theorem, which is one of the principal results of the paper, under an additional assumption 

on the smallness of the capacities of F*^^-* . 

Theorem 1. Let the parameters e and d, introduced in (5), satisfy the inequality 

e<c (53) 

where c is a sufficiently small absolute constant. 

Then the matrix I + SD, defined according to (28), (29), is invertible, and the solution m(x) to the 
boundary value problem (8)-(9) is defined by the asymptotic formula 

N 

u(x)=t;/(x) + ^Q(p(^')(x)-47rcap(J^(^')) if(x,0(^'))) + i?(x), (54) 

where the column vector C = (Ci, . . . , Cjv)^ is given by (30) and the remainder -R(x) is a function harmonic 
in flN, which satisfies the estimate 

|ii(x)| <c{6\\Wvfh^^^^+e'd-y'\\vfh^^^)]. (55) 



Proof: 

The harmonicity of R follows directly from (54). 
If X G dCl, then 

i?(x)= -4.^Qcap(F«)(-^^-^-if(x,0«)) 

JV 

+ ^\Cj\0{e cap(i^W)|x-0(^)r2). 
Since G(x, O^^)) = on dn, and P(j) satisfies (15) we deduce 

N 

i?(x) = ^0(e cap(p(^))|Q||x- 0«|-2), (56) 

where |x — O^-'^ \>Cd, and C is a sufficiently large constant. 
If X e a(M3 \FW) then 

N 

i?(x) = -^;/(0«) + 0(£||V^;/|U„(^)) +47rX]C,cap(F«)(iJ(OW,0«) + 0(£)) 

i=i 

^ r cap(F(^)) £ cap(FO)) -( 
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Noting that (30) can be written as the algebraic system 

N 

-fl"(0('=\ 0(J'')) + t;/(OW) = 0, (58) 

which, along with (57) and the obvious inequaUty cap(i^*^-')) < e, impUcs 

i?(x) = 0(£||Vt;/|U^(„)) +47rCfecap(F«)iJ(0W,0W) 

+ 5:0(.cap(F«)|C,|)+ Yl Q( |o(r-Q(i)|2 l^^-|)- (59) 

j=l l<j<N, jjik I I 



It suffices to estimate the sums 

E 

and 



£ cap(f (J))|Q| 

lo(fe) - owp 



^ lx-0«P ' 

l<j<N ' ' 

When £ < c d^/^ we refer to Lemma 1, and using the inequality (32) we derive 

1/2 / ^ 1/2 



^ £cap(FO))|Q| / ^ ^^cap(fO)) \ .^0)^^2y 
3¥'k,l<j<N ' ' \j¥'k,i<j<N ' ' / \i<j<iv / 

< const 1^ ^ cap(i^(^))(t;/(0(^)))2j ^ ^niax^£2d-3cap(F(j)) 



i<j< 

£2 

< const (60) 
Similarly, when x G 90 we deduce 

£ cap(J^(^))|C,-| £2 „ „ 

E , _o0)|2'' < const (61) 

l<j<N ' ' 

Combining (60), (61), (56) and (59) we complete the proof by referring to the classical maximum principle 

for harmonic functions. □ 

Under the stronger constraint (43) on e and d, Lemma 2 and representations (56), (57) lead to the 
following 

Theorem 2. // the inequality (53) is replaced by (43), then the remainder term from (54) satisfies the 
estimate 

|i?(x)| <c{e\\Vvfh^^^)+s'd-^Vfh^^^)}. (62) 
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7 The energy estimate 

Under the constraint (43) on e and d, which is stronger than (53), we derive the energy estimate for the 
remainder R. This result is important, since it allows for the generalization to general elliptic systems, and 
in particular to elasticity where the classical maximum principle cannot be applied. 
Theorem 3. Let the parameters e and d, introduced in (5), satisfy the inequality 



where c is a sufficiently small absolute constant. Then the remainder R in (54) satisfies the estimate 



£ < c d^, 



(63) 



||Vii|U,(o,) < Const 



(64) 



Proof. 



For every k = 1,...,N, we introduce the function 




(65) 



JV 




AnCk cap(FW)i?(0('=),0('=)) 



where the coefficients Cj satisfy the system (26). 




i?(x) + *fc(x) = -z;/(OW)-Cfc 
V- ^ f cap(F(.')) 




which together with (26) implies 




i?(x) + *fe(x) = 
outside of a set with zero capacity). 



The function ^q, defined by 




N 



(66) 



satisfies 



i?(x) + *o(x) = 

quasi-every where on dft, which follows from (54) and (9), (11). 
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Wc set B^''' = {x : |x — O'^'^^l < p}, and define the capacitary potential of F^''^ relative to -B^^^, that is 
a unique variational solution of the Dirichlet problem 

APfe(x) =0, xeS(5^\i^W, (67) 

Pfc(x) =1, xea(M3\F('=)), (68) 
Pfc(x) =0, |x-OW|=rf/4. (69) 

Also, let a surface Sd be a smooth perturbation of dQ, such that 

SdCfl and rf/4 < dist(5d,x) < d/2 for all x G Ofl. 

In turn, the set of all points placed between the surfaces and Sd is denoted by 11^, and the function Pq 
is defined as a unique variational solution of the Dirichlet problem 

APo(x) =0, xGHd, (70) 
Po(x) =1, xeaf2, (71) 
Po(x) =0, xG5d. (72) 

We note that 

N 

P(x)+^Pfe(x)vl/fc(x) (73) 
fe=o 

vanishes quasi-everywhere on dflN and that the Dirichlet integral of (73) over fijv is finite. Therefore, by 
harmonicity of R 

I VP(x)- V(p(x)+ A(x)*fc(x))dx = 0. 

"'^JV 0<k<N 

Hence 

l|Vi?||L(a.) < l|VP||L.(n„) ||V Pk'i^kh.in.), 

0<k<N 

which is equivalent to the estimate 

E II V(Pfc*fe) f + II V(Po*o) |lL(n.) • (74) 

fe=i I 

In the remaining part of the proof, we obtain an upper estimate for the right-hand side in (74). 
The inequality (74) and the definition of lead to 

II VP ||i,(a,) <2(/CW+/C(2)+z:(i)+£(2)+x(i)+A^(2)+^+q), 

where 

x;w =Ellv(^fc(^/(-)-^/(o('=))))||^^(^(.)), (75) 
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k=l l<j<N, jjik ' ' 

AT N 

= (47r)2^ II ^C, cap(F(^)) V (^Pk{H{-,0^^'>) 
fe=i j=i 

-/f(OW,OW))) ' (77) 

N 2 2 

AT =(47rf ^ICfcl^ (cap(F('=))) (^(0^,0^)) || VPfe||^^(^(.K, (78) 



fe=i 



Q =11 E '^.^('^.(^'"0-g^))llL,n„. m 

' are defined by repl 
We start with the sum K.^^^ . Clearly, 



and /C(2), £(2), 7W(2) are defined by replacing B^^^ in the definitions of /C^^), £''^\ M^^'^ by B^/^ \ B^^^ 



AT 



<t?l|V«/||L(,) ^ / {|VP,(x)|2 |x-OW|2 + (p,(x))'}dx 

<^?I|V«/||L(^) f^e^cap(pW) (80) 
fe=i 

and hence 

< C£3d-^||Vz;/||i^(„). (81) 
Furthermore, by Green's formula and by (10) we deduce 

N „ 

fc=l -^^d/A^Se 

+2VPfe(x) • Vw/(x)|dx 
-E P.(x)(.,(x) -./(0('=))){p.(x)|g(x) 

+ (^;/(x)-^^/(0«))|^(x)}d5 (82) 
By the mean value theorem for harmonic functions and the inequality Pfe(x) < P('^'(x), we have 

where B = {y : |y — x| < |x — 0('^^|/4}. Making use of the asymptotics (15) far from O^*^^ we deduce 
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Now we turn to the estimate of (82). The volume integral in the right-hand side of (82) does not exceed 
C i:ilV.,|U.<„,cap(F«) / {^^l/MI 

fe=l ■'^d/4\^3e I I 

cap(F('=)) „^ „ ^ , 



< 



<Cs^rf-^ll/llL(n.)- 
By < Pfe(x) < 1 and (83), the surface integral in (82) is dominated by 

Combining (84) and (85) we arrive at the estimate 



-3^-3|| f\\2 

llL„(nN)- 



Let us estimate C^^'> (see (76)). Obviously, 



fe=l \ l<j<N, jjik ' 



. l<j<JV, j^k 

Furthermore, when j ^ k we have 



< 



which does not exceed 



cap(F(j)) 

10^=) - oo) 



£ cap(ii^(j))(cap(F('=)))V2 ^ £3/2(.ap(^(i)) 



low - OW|2 



|0(fe) - O0')|2 



Hence, using Lemma 1 we deduce 



fe=l \ l<i<Af, j^k ' ' / 

N 



fc=l l<j<N, j^k 



l<J<iV, jy/c 



and therefore 



(84) 



(85) 



(86) 



(87) 
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Similar steps can be followed to estimate Q in (79). We have 



< 



cap(i^(j')) 
|x-O0)| 



which does not exceed 



i2(n<i) 
i2(n<i) 



P„(vp(^)+cap(F(i))^^) 



i2(nd) 



Ce cap(F(^))|x - O^^'^T^, x e Hd. 
As above, we use Lemma 1 to deduce 

2 



\l<i<Ar ' 



1 



l<j<Af l<j<N 

Next, we estimate C^'^\ Integration by parts gives 



-O0')|4 



<C£^d-^||t;/||i^(^) for xGHrf. 



|OW - O0)|. 

j(i^a)) 



(88) 



^^«(-'^Hx)-^i^){P.(x)^P-(x) 

-(--(x)-^i^)l^(x)}... 
When j 7^ fc and k, the volume integral in the right-hand side of (89) is estimated as follows 

P,(x)(pW(x) - j^^^ (VP.(X) . VPW(x))rfx 



(89) 



< C 



< C 



cap(P('=)) cap(P(^)) 



-O0)|2 A 



cap(P(*^)) cap(p(™)) 



bW |x - OC^) |2 |x - 0('") |- 



dx 



(90) 



< C- 



cap(P(j)) cap(p('")) 



(cap(p('=)))2cap(p(^))cap(p(™))d 



In turn, when j ^ k and m ^ k the modulus of the surface integral in the right-hand side of (89) does 
not exceed 



C 



e cap(F(-'')) 



cap(p(™)) ^ ircap(p('")) SPfe 

ail (x) 



>dS 
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^ e cap(i^O-)) cap(F("')) \ ^ ^ f cap(fW) \ 
- |OW-O0)|2|OW-O(™)|2|^ iasw Ix-OWp"*^] 

^ £3 cap(FO-)) cap(f (-)) 



We have 



i<m,j<N l<k<N, kjim,kjij '^d/4\'^3e ' ' 

and by (89) , (90) and (91) 

r(2)<rp2, V- ir \\r \ V cap(i^(j))cap(F("')) 

= ^^5 E l^™||^?.|cap(i^(^))cap(F(")) ^ |o(^) - QOOHoC^) _ QMp - 

i<m,j<N i<k<N, k=^m,k^j ' ' ' ' 

and therefore 



r(2)<r^" V- |C„.||Q|cap(F(.'-))cap(f('")) 

l<rn,j<N ' ' 



Let us introduce a piece-wise constant function 



C(x) = ( |Cm|(cap(f W))V2, when x e 

1 0, otherwise. 



(m) 

d/4 ' 



Then the inequahty (92) leads to 

(2)^ £3 (|C^|(cap(fM))V2) (|Q|(cap(f-(^)))V2)d^ 

l<m,j<N ' ' 

where the constant C depends on lo, and using Lemma 1 we deduce 

^^'^^^i E C'|cap(F«)d3<c (93) 

To evaluate M^^^ + M^"^^ we apply the result of Lemma 2 and use the same algorithm as for K.^^^ and 
/C(2) to deduce 

A^(i) +^(2) < C\\vf\\l^^^^ed-'{e'd-'+e'd-') < Ce'd-'\\vf\\l^^^y (94) 

Similarly, applying Lemma 2, we derive the estimate for the term Af 

^f<Ce^r'\\vf\\l^^^y (95) 
The proof is completed by the reference to (80), (86), (87), (88), (93) , (95). □ 
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8 Meso scale approximation of Green's function in Q^v 

Let Gjv(x, y) be Green's function of the Dirichlet problem for the operator —A in Ojv- In this section, 
wc derive the asymptotic approximation of GAr(x, y) and estimate the remainder term. In the asymptotic 
algoritlim, we will refer to the algebraic system similar to that of Section 5. We need here Green's functions 
g^^^ (x, y) of the Dirichlet problem for the operator —A in \ F^^^ , j = 1, . . . , N. The notation h^^^ will be 
used for the regular part of ^(•'^ that is 

/i(j')(x,y) = (47r|x-y|)-i-5W(x,y), x, y e \ F^^'^ . (96) 

According to Lemma 2 of [2], the functions h^^^ allow for the following estimate: 



x-O0)P' 



for all y e \ F^) and |x - O^J' | > 2e. 

The principal result of this section is 

Theorem 4. Let the small parameters e and d, introduced in (5), satisfy the inequality e < c , where 
c is a sufficiently small absolute constant. Then 

Gn (x, y) = G(x, y) - 5^ I /i^) (x, y) (98) 

-P(^'Hy)^(x,0(^')) -P(^)(x)iJ(0(^'\y) +47r cap(p(^'))ir(x, 0(^))ir(0(^'), y) 

N \ 
+if(0(^),0(^)) T(^)(x)rW(y) - ^Ci,r«(x)T(^)(y) \ +7^(x,y), 

i=i J 

where 

T(^\y) = P(^')(y) - 47r cap(F(^))iJ(0(^\y), (99) 
with the capacitary potentials P^^^ and the regular part H of Green's function Goffl being the same as in 

C = (I + SD)-^S, (100) 



Section 3. The matrix C = (Cij)f.-^i is defined by 



where S and D are the same as in (28), (29). The remainder 'TZ{x, y) is a harmonic function, both in x and 
y, and satisfies the estimate 

|7^(x,y)| < const ed"^ (101) 

uniformly with respect to x and y in flN . 

Prior to the proof of the theorem, we formulate an auxiliary result. 

Lemma 3. Let the small parameters e and d, defined in (5), obey the inequality (43). Then the matrix 
C in (100) satisfies the estimate 

||c||r«^m^ < cc^-^ (102) 

where c is an absolute constant. 
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Proof. First, we note that 

l|i^llR-^'^R« < const ||S||ra'^r«, (103) 

which follows from Lemma 2, where V/ should be replaced by the columns of the matrix S. 
Additionally, 

||S||Riv^Riv < const d~^. (104) 
To veriiy this estimate we introduce a vector ^ = {^j)jLi, \\^\\ = 1, and a function ^(x) defined in Q by 

C(x) = I = ^ I" \<d},j = l,...,N, 

' '0 otherwise. ^ ' 



Then 



(SCC) < const / / ^(X)eCy)'^XdY ^ ^^^^^ j |^(x)|2rfx < const d-^ 

Jc^Jc^ 47r|X-Y| 

which yields (104). Then (103) together with (104) lead to (102). □ 

Proof of Theorem 4- The harmonicity of TZ follows directly from (98). 
Let us estimate the boundary values of TZ on O^n- 

If X G dQ and y & ^n, then according to the definitions of Section 3.3 for Green's function of O and its 
regular part the remainder term TZ in (98) takes the form 



+H(0«),0«>) TU)(y)(p«)(^) _ 



N 

-^c.,r(^)(y)(p«(x) 



cap 



(i^W) 



lx-OW| 



Taking into account the estimate (97) for /i^-'^ together with the asymptotic representation (15) of P^^^ we 
obtain 

N r>(i)f \ N N 2 

«(-.^) - Eo(3^) +i:E^«^«'(y)o(^3WF) f'""' 

j=l I I j=l i=l I I 

for all X G dfl. 

Here 7?.(x, y) is harmonic as a function of y. Next, we estimate (106) for y G OCIn- 
If y,x G dn then (106), (99) and (15) lead to 

N 2 N N 4 

Using (102) we can estimate the double sum from (107). For a fixed x G dfl, let us introduce a vector 

V|x-OW|2y'i=i' 
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and a function V{X.) defined in fl by 

V{X) 



Vj when X G B(^), 
otherwise, 



where the balls B^^^ are the same as in (105). It follows from Lemma 3 that the double sum in (107) does 
not exceed 

^""'^ / (nx))^dx < 



The above estimate together with (107) imply 

7^(x,y) = 0{£^d-^\ \ogd\ + e^d"^) when x,y e dO.. (108) 
Now, we estimate (106) for y e d{R^ \ F^™'). In this case we have 

_2 



^("'^)= ^( |x-O0)p )+ ^ Q( |y_oo-)yx-oo)p ) 

^ 2 f 

+ ^^( |x-OWp )^-(^-^" cap(FW)if(OW,y)) 



+47r J2 cap(f ^^'))G(0<^\y) +Q( |y _ Q(^.)p )) }• (109) 



We also note that according to (100) the coefficients Cij satisfy the system of algebraic equations 
(5i„)G(0('"),0«)-Ci„-47r ^ cap(F«) G(0("), 0«) = 0, m, i = 1, . . . , TV. (110) 

l<j<N, j=im 



Hence, in the above formula (109) the expression in curly brackets can be written as 

.2 

r 

l<j<N,j^m 

and then formulae (109) and (111) imply 



Cim + 0{\C,m\e) + ^n Yl C,,(cap(FW)G(OW,y)+0(— ^— )) (111) 

= (i-5,„)G(o('"),o«) + 0(£d-^)+ ^ qMt-^^), yG9(K^\i=^('")), 



7e(x,y)=o(erf-^+.^|logdM-^+.^d-^)+ ^ Y: ^-^•Q( |v-O0)|C-O(-)p )' 

l<i<N l<j<N,j^m ' ' ' 

where the estimate of the double sum is similar to (107). Thus, we obtain 

7e(x,y) =o(£d-2 + e2|log(i|rf-3 + e3d-^ + e4(/-7^ =0{ed-^), (112) 

for all X e 51^ and y e d{R^ \ F('")), m = 1, . . . , iV. 
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Using the estimates (108) and (112) and applying the maximum principle for harmonic functions we 
deduce that 

7e(x,y) =0(£d-2), (113) 

for all X G dO, and y G 

In turn, when x <E 9(R^ \ F^'''' ), the formula (98) and the definition (96) of h^^^ lead to the expression for 
the remainder term on the boundary of the inclusion 



7^(x,y)= i?(x,y)-ir(0«,y)+ ^ (/iW(x,y) - P«(x)i?(0«,y)) 

N / N \ 

+ ^T(^)(y) j ir(o(^), 0(^)) t(J)(x) - if(x, 0(J)) - ^CyTW(x) . (114) 
j=i \ i=i / 

Using the formulae (15) and (97) for P^^^ and h^^^ together with the definition (99) of T^^^ and the definition 
of Section 3.3 of the regular part of Green's function of fl we deduce that 



(115) 



47r|x ■ 
and 

i/(x,y)=iJ(0W,y) + 0(£), (116) 
for X € a(R3 \ F^^)) and y € 17^. The representations (99) together with (114)-(116) imply 



^(^'^)= ^ 4.|x-O0-)| +^(p^ 

JV / 

_ ^T(^)(y) iJ(oW, 0(^)) - iJ(0(^\ 0(^))T(^')(x) 



+ f]C,,TW(x)] +f;0(£|TW(y)|). (117) 

Bearing in mind the asymptotic formula (15) for the capacitary potentials and the definition (99) we deduce 
that for X e ^(R^ \i^(fc)) 

2 

= 47r cap(F«)G(x, O^^)) + o(^--^^) , A:. 
Thus, (117) can be rearranged in the form 

7e(x, y) = ^ r(-') (y) 1 0(0^=) , 0(^") ) - C^j 

-47r ^ Q,- cap(FW) G(OW,0«)l +7^W(x,y), (119) 

l<i<N, i^k J 
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where 



l<j<N, jjik ' ' 

^ ^ £cap(F«) 



low -0«|2 

l<i<JV, iT^fc I I 



+1: E 



e cap(i^(*)) 

It follows from (100) that the coefficients satisfy the system of algebraic equations 



(1 - <5fc,)G(0(^), 0(^)) - Ckj - 47r Yl cap(F«) G(OW, o«) =0, = 1, . . . 

l<i<iV, ijtk 

and hence using (119)-(121), we arrive at 

7e(x,y)=7eW(x,y) 

for all X G \ F'^''^) and y G f^^. 

Let us consider the case when y G d{M.^ \ F^™)). Then 

TW)(y) = 47r cap(FW) G(0(^-),y) + o(i^^P^), j ^ m, 

and 

T(™)(y) = 1-477 cap(F(™)) ff(0('"\y). 
The double sum in (120) can be rearranged according to (110) 



,|0(fc) - OW|2. 



~t~ ^(l^zml^) 



= y ecap(FW) x f 

l<i<N, ijik ' ' l 

+47r ^ C., ( cap(f )G(0(^) ,y)+0{ / ))] 

l<i<N, ijik,ijtm ' ' 

= 0{e^\logd\d-^ + e^d-^ + e^d""^), for x e ^(IR^ \ F^^)), y g ^(IR^ \ i^W), 
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where the estimate of the last double sum in (123) is similar to (107). Combining (120), (122) and (123), 
we deduce that 7^(x,y) = 0{ed-^) for x G ^(R^ \ fW), y e d{m?' \ F^"^^), m,k = 1,...,N. Using the 
symmetry of H{x,y) together with (112) we also obtain that 7^(x,y) = 0{sd-'^) for x G d{R^ \ F'^''^), k = 
1, . . . ,N, y G dil. Applying the maximum principle for harmonic functions we get 

7^(x, y) = 0{ed-'^) for x G d{m^ \ F^'^^), k = l,...,N, y G f2jv. (124) 

Finally, formulae (113) and (124) imply that TZ{x,y) = 0{ed~'^) for x G 90jv and y G Ojv, and then 
applying the maximum principle for harmonic functions we complete the proof. □ 
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